Janssen effect and the stability of quasi 2-D sandpiles 
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We present the results of three dimensional molecular dynamics study of global normal stresses in 
quasi two dimensional sandpiles formed by pouring mono dispersed cohesionless spherical grains into 
a vertical granular Hele-Shaw cell. We observe Janssen effect which is the phenomenon of pressure 
saturation at the bottom of the container. Simulation of cells with different thicknesses shows that 
the Janssen coefficient k is a function of the cell thickness. Dependence of global normal stresses as 
well as K, on the friction coefficients between the grains (/^p) and with walls ( ^^) are also studied. 
The results show that in the range of our simulations n usually increases with wall-grain friction 
coefficient. Meanwhile by increasing /ip while the other system parameters are fixed, we witness a 
gradual increase in k to a parameter dependent maximal value. 

PACS numbers: 



I. INTRODUCTION 

Pouring a large amount of granular material into a 
vertical and static cell produces two different, interesting 
structures depending on the geometrical boundary con- 
ditions imposed on the system. In a silo when all the side 
walls are closed, injection of granular material produces 
a static granular pack which is a vertical column with 
a flat and horizontal free surface. Granular packs have 
been the subject of studies for long timoir— because, in 
addition to the practical applications like corn storages in 
silos, they serve as a useful model for a variety of physical 
systems. On the other hand, inserting grains into a cell 
with at least one open boundary, generates a stable heap 
with finite slope. This non-zero angle to the horizontal, 
the angle of repose 0, is one of the most important charac- 
teristic of the granular material at macroscopic scales^"—. 

A remarkable feature of deep granular columns is the 
saturation of vertical stress at the bottom of the silo 
known as the Janssen effect, named after German en- 
gineer H.A. Jansseni. Janssen was one of the first who 
studied systematically the effect of pressure screening at 
the bottom plate of silos and succeeded to present a theo- 
retical framework which accounts for this peculiar behav- 
ior of granular matter. In Janssen's model the granular 
material is treated as a continuous medium in which a 
fraction k of the vertical stress cr^, converts to the hor- 
izontal stress. With the further assumption that all the 
frictional forces exerted by the walls are at the Coulomb 
failure criterion, he derived the following expression for 
the value of Uzz at the depth z of a granular pack in a 
vertical container: 



,(z) = pgl{l 
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where p is the mass density of the granular material, 
g is the gravitational acceleration, and I is the decay 
length. The decay length (and subsequently the satu- 
rated stress cr™"^ = pgl) is related to the cell geometry 
and wall-grain interaction as well as the grain-grain in- 
teractions. Janssen's theoretical analysis leads to a value 
of I = 77^^ for decay length in which Ai, and U are 
the area and circumference of the cell profile respectively 
and pw is the friction coefRcient between the walls and 
grains (sometimes and especially in old literature, the 
value oi K = nfj^w has been referred to as Janssen's coef- 
ficient). Recent experimental and computational studies 
revealed that the above formula needs to be modified, 
and as such, the classical Janssen analysis and specif- 
ically the assumption of Coulomb failure is not always 
satisfied. A number of modified Janssen formulas (which 
usually mandate a region of perfect hydrostaticity next 
to the surface) have been suggested which fit the data 
obtained from real experiments and sophisticated siniu- 
lations^iii^r— better than the standard Janssen's formula. 
However, the concept of pressure saturation at the bot- 
tom of the container is still valid. 

The Janssen effect originates from the ability of gran- 
ular materials to support shear stress. The frictional 
forces exerted by container's lateral walls carry part of 
the grains weight. These very frictional forces are also 
responsible for the formation and stability of a sandpile 
in a container with an open boundary. At the begin- 
ning, when the injected grains fall on the bottom of the 
cell, the frictional forces between the grains and the cell 
walls dissipate the major part of the grains kinetic en- 
ergy, helping the rapid formation of a sandpile. Mean- 
while, the presence of lateral walls modifies the angle of 
repose mostly because of the arching phenomena between 
front and rear walls^. In fact, in the case of a thin con- 
tainer (usually called the granular Hele-Shaw cell) many 
studies have demonstrated that just like the saturated 
normal stress, the angle of repose 9 is dependent on both 
cell geometry (here the cell thickness w) and the wall- 



grain friction cocfHcient^rii. Courrech et al — proposed 
a simple model to incorporate the lateral wall effect on 
piling via the Janssen effect. This model introduces a 
new characteristic length defined as Br = 2Kfiwhfreeze, 
where h freeze is the flowing layer height when an induced 
surface avalanche is about to stop. The predicted growth 
of 9 in Courrech ct al's model: 



Tr\ "^"^"(^^^ "' ) 

COS[0^) Br 

is slower than the empirical exponential law^: 



(2) 



(3) 



which has been traditionally used to fit their data. Here, 
a is a constant which depends on grain properties and 
A is a characteristic length representing the scale over 
which the walls affect the piling of grains. 

In this paper we are interested in numerical evaluation 
of the applicability of the Janssen assumptions in the 
case of quasi 2-D sandpiles formed by pouring granular 
material into thin vertical Hele-Shaw cells. First we show 
that, similar to granular packs in silos, the global pres- 
sure at the cell's bottom, at, saturates to a limiting value 
while the mass is increasing. We then investigate the ef- 
fect of wall-grain and grain-grain friction coefficients on 
(Tb and also global horizontal stress ay at the front and 
rear walls. An important issue is the dependence of the 
Janssen coefficient n on system parameters. We examine 
the effect of the friction coefficients Hp and /i-u,, as well 
as cell thickness w, on the value of the Janssen coeffi- 
cient and observe interesting and sometimes unexpected 
results. 
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FIG. 1: A schematic view of the simulation box 



velocities v,; , Vj and angular velocities uji , cuj the force on 
grain i is computed as follows. The normal compression 
6ij is 



Si, J = d-Ti, 



(6) 



The relative normal velocity v„. . and relative tangential 
velocity vt . . are given by 



{^r^,j.n^,J)n.i, 



(7) 



II. MODEL 



vt, J. = Vij - v„,^ - -(wj + t^j)rj.j 



(8) 



Our simulations are based on a Distinct Element 
Method (DEM) scheme, originally proposed by Silbert 
et ali^. For the sake of clarity, we review the most im- 
portant aspects of it here. According to this model, the 
spheres interact on contact through a linear spring dash- 
pot interaction in the normal and tangential directions 
to their lines of centers. In a gravitational field 5, the 
translational and rotational motions of grain i in a sys- 
tem at time i, caused by its interactions with neighboring 
grains or walls, can be described by Newton's second law, 
in terms of the total force and torque with the following 
equations: 



Ff* = m,g + ^(F„,.,+F,„ 



(4) 



where n. 
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with r.ij — |rjj| and v^ j — -Vi — Wj. The 



rate of change of elastic tangential displacement Uj^ . , set 
to zero at the initial part of contact, is given by 



dxif 



dt 



Vt. 



(ut.,,-v,,j)r, 



(9) 



The second term in Eq.(9) arises from the rigid body 
rotation around the contact point and insures that Ut . . 
always lies in the local tangent plan of contact. In Eqs. 
(1) and (2), the normal and tangential forces acting on 
grain i are given by 



^"i,i = fi^)ik„Si,jnij ~ 7n"^e//V„, J (10) 



T; = — 
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Where tti, and r. 
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Yi — Yj are respectively, the mass 
of grain i and the relative distance between grains i and 
j. For two contacting grains i, j at positions r^, Yj with 



and 



F*.. = f{^){ktnu.,-ltmeff^u,,) (11) 

Where fc„.t and 7„.t are elastic and viscoelastic constants 
respectively and meff = mimj/(rni + mj) is the effective 



mass of spheres with masses mi and rrij . The correspond- 
ing contact force on grain j is simply given by Newtons 
third law, i.e., Fij ~ —Fj^i . For spheres of equal mass 
m, as is the case in our system, rrieff = m/2; f{x = 1) 
for the linear spring dashpot (Hookian) model with vis- 
coelastic damping between spheres^. 

Static friction is implemented by keeping track of the 
elastic shear displacement throughout the lifetime of a 
contact. The static yield criterion, characterized by a 
local grain friction coefficient ^, is modelled by trun- 
cating the magnitude of u^^ - as necessary to satisfy 
|F(.^. I < |/iF„; J. Thus the contact surfaces are treated 
as sticking when |Ft.^. | < |^F„; J, and as slipping when 
the yield criterion is satisfied^s^S. 

In the following we present the results of our extensive 
molecular dynamics (MD) simulations^! in three dimen- 
sions on model systems of N mono-disperse, cohesionless 
and inelastic spheres of diameter d and mass m. The 
system is constrained by a rectangular box with fixed 
rough walls and free top surface, as in figl. A simulation 
was started with the random generation of spheres with- 
out overlaps from top and left corner of container, fol- 
lowed by a gravitational settling process to form a stable 
heap. The results are given in non-dimensional quantities 
by defining the following normalization parameters: dis- 
tance, time, velocity, forces, elastic constants, and stress 
are, respectively measured in units of d, to = \/d/g, 
Vo = \/dg, Fq = mg, fco = "^ff/c^i a-nd a = mg/d'^. All 
data was taken after the system had reached the steady 
state. Because of the complexity of the model, there are 
a wide range of parameters that affect the results of com- 
putation. We usually investigate the effect of a single pa- 
rameter varying in a certain range while other variables 
are fixed to their base values as listed in Table HI 

All the cases were simulated in three dimensions us- 
ing a molecular dynamics code for granular materials 
LAMMPS^i2^. The equations of motion for the transla- 
tional and rotational degrees of freedom are integrated 
with either a third order Gear predictor-corrector or 
velocity- Verlet scheme^^. 



TABLE I: Basic computational parameters 



Parameters 



values 



Maximum Number of grains (A'^) 
grain-grain friction coef . (pp) 
wall-grain friction coef. (p-w) 
grain normal stiffness coef. (fc„) 
grain tangential stiffness coef. (kt) 
grain normal damping coef. (7^) 
grain tangential damping coef. (74) 
wall normal stiffness coef. (fc„) 
wall tangential stiffness coef. (kt) 
wall normal damping coef. (7^) 
wall normal damping coef. (71) 
Time step increment 



40, 000 

0.5 

0.5 

2 X 10=* (fco) 

2/7fc„ 

50/(io) 

50/(io)) 

2 X 10=* (fco) 

2/7fc„ 

50/(to) 

50/(io) 

2 X 10"=* 




FIG. 2: The geometry of the pile and the wall normal stresses. 



III. SATURATION OF WALL NORMAL 
STRESSES WITH MASS 

We start from an empty vertical cell composed of two 
parallel plates (front and rear walls) separated by an 
spacer with an adjustable width w. A granular heap is 
formed in the cell by pouring the grains on the bottom 
plate. These grains are released from a small box located 
on the top-left of the cell (see fig.l and fig. 2) where they 
have been randomly located in the first place. The num- 
ber of grains filling the small box, depends on the ma- 
terial density and the size of box. The rate of material 
insertion R can be varied by changing the box size. The 
process repeated successively. A pile with well defined 
shape starts to form as the number of the added grains 
grows. The angle of repose can then be determined from 



the surface profile of the heap. The angle of repose start 
from zero and reaches very soon a stable value. Previ- 
ous studies have demonstrated that the angle remains 
constant with further insertion of grainsi^. 

The growth of the height of the pile H and length L 
(measured from the pile profile) with grain number N in 
a typical cell with size w = 6 is shown in fig. 3. All the 
other parameters are as those listed in Table HI Clearly, 
the total pile mass M, is proportional to the number of 
grains in the cell. It can be seen from fig. 3 that after a 
short initial time, both L and H grow linearly with y/N. 
This is the expected behavior of a pile with constant 
density p and fixed angle of repose 0. 

The total normal forces exerted by the front (or rear) 
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FIG. 3: Variation of the Length L (triangles) and the height 
H (circles) of the sandpile with the square root of grain num- 
ber A*' in a cell of thickness w = 6 



FIG. 5: Variation of the front (or rear) wall normal stress (Ji, 
(triangles), bottom pressure at (squares), and normal stress 
on the spacer CTs (circles) with grain number A'^ in a cell of 
thickness lo = 6. 
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FIG. 4: Variation of front or rear wall normal force, Fv (trian- 
gles), apparent weight. Ft (squares), and normal force on the 
spacer, Fs (circles) with grain number Ai" in a cell of thickness 
ui = 6. 



FIG. 6: Estimation of Janssen coefficient k — av/cTb from 
the data presented in fig. 5 in a cell of thickness w — 6 with 
iip = fx-w = 0.5. 



wall, F^, the spacer, Fg^ and also the apparent weight 
of the grain pile, i.e. the normal force exerted by the 
bottom plate, Fb, are all increasing functions of mass as 
depicted in fig. 4. Further investigations of these curves 
suggest that while F„ grows linearly with N , both F^ and 
Fs increase linearly with ^/N. Given the fact that except 
for the thickness w the other pile linear sizes are growing 
as V^, one can deduce that all the global normal stresses 
(Jb = jj^ (bottom pressure), (Js ~ -^ and a^ — ^j^ (hor- 
izontal stress) should approach constant values after the 



transition period, that is when a well-defined heap has 
been formed (sec fig. 5). From the data presented in fig. 5 
one can see that compared to cr„ , the value of as is larger 
at the beginning of the pile formation process. In other 
words, in small heaps the horizontal stress in the pile is 
not isotropic. However, as anticipated for very large val- 
ues of iV, Os eventually becomes almost equivalent to cr„, 
as in to the case of 3-D granular columns. 

From the obtained values for normal stresses, we have 
calculated k, the fraction of bottom pressure Ub which is 
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FIG. 7: Variation of bottom pressure Oh with grain number 
Af in a cell of thickness u; = 6 for different values of ^Lw'- 
11^ = 0.1 (□); Mu. = 0.3 (O); ^» = 0.5 (A); ^^ = 0.7 (V); 
M» = 0.9 (+). 



FIG. 8: Variation of horizontal cr„ with grain number N in 
a cell of thickness w = 6 for different values of /iw'- fJ.w =0.1 
(□); ^» = 0.3 (O); /i„ = 0.5 (A); ^» = 0.7 (V); /i„ = 0.9 
(+)• 



taken up by the front and rear walls, ay, as a function 
of heap mass. The results presented in fig. 6 demonstrate 
that , starting from smaller values, the Janssen coeffi- 
cient K reaches very quickly to a maximal final value 
during the filling process. This is in accordance with 
the Janssen theory which assumes a constant value for n 
independent of the mass of the injected material. Using 
the evaluated values of k in the flat region of the curve, 
we have estimated k = 0.51 ± 0.01 for this very case. 



IV. EFFECT OF FRICTION COEFFICIENTS 
ON WALL NORMAL STRESSES 



The difference between real and apparent weight in 
granular packs and confined sandpiles stems from the 
ability of the medium to redirect part of the normal stress 
experienced at the bottom of the container to the vertical 
walls surrounded it. For this we need a stable network of 
grains in contact with each other which spreads the forces 
through at the medium. Obviously, the frictional forces 
plays the central role here, as they sustain the stability 
of the contact networks and heap morphology. 

Generally, the contribution of side walls, i.e. the shear 
stress they impose on the confined granular material de- 
pends on the strength of wall-grain frictional forces which 
is determined by fi^ , the wall-grain coefficient of friction. 
To examine the behavior of the system when fi^ changes, 
we performed a series of simulations on different heaps 
all formed in a cell of thickness w = 6 but with different 
wall friction coefficient fiw All the other parameters in- 
cluding grain-grain friction coefficient have been fixed to 
the basic values listed in Tabled) From the data we have 




FIG. 9: Variation of bottom pressure ai, versus ^^ for A'^ 
40000 in a cell of thickness w = 6, with fip = 0.5. 



calculated the variation of normal stress ab, and the be- 
havior of the horizontal stress a^, , with the grain number 
N. For more clarity we have also sketched in separate 
sheets the values of ai, and ay as a function of fi^, for 
the maximum grain number wc used in our simulations, 
N = 40000 (sec fig.9 and fig. 10). Inspection of these 
curves shows that the wall normal stresses almost cease 
to change with ^^j when it passes the value oi ^w = 0.5, 
reach to their final minimal values. In other words, for a 
granular material with specified grain parameters, there 
is an upper limit to the contribution of walls in carrying 
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FIG. 10: Variation of horizontal stress (t„ versus /i„ for N 
40000, in a cell of thickness w = 6, with fip = 0.5 
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FIG. 11: The dependence of k, on /ii„ in a cell of thickness 
w = 6, with /ip — 0.5. 



the weight of material. The case of iip — 0.1 is an ex- 
ception. The frictional forces between the container (and 
specially the friction applied by the bottom) are too small 
and as such the formation of a heap with steady angle of 
repose needs more times (grains). 

From the above data wc are also able to make an al- 
most precise estimate of the dependence of k on /i„ by 
calculating the value of a-ujah for different values of [i^ 
(fig. 11). Keeping track of variation of k during the filling 
process we find out that the fraction of bottom pressure 
converted to horizontal stress is almost constant, and the 
variation of k with mass always follows a trend like that 



of fig. 6. Therefore, the estimated values for Janssen co- 
efficient, presented in fig. 11 are precise even for the case 
of /i^ = 0.1. This figure also suggests that k grows with 
\Xw but the growth slows down at larger values of [i^ ■ 

It is believed that one of the most effective medium pa- 
rameters at the microscopic level is the grain-grain fric- 
tion coefficient yi,^^. Therefore, we performed a series of 
simulations for some selected values of /ip while the other 
parameters have been kept constant. The obtained val- 
ues for wall normal stresses g\, and a^ in a pile consisting 
of at most TV = 40000 grains have been presented in fig. 12 
and fig. 13 respectively. Comparison of fig. 12 with fig. 9 
shows that increasing each of the friction coefficients (^p 
or ^uj) decreases Oh to a parameter dependent asymp- 
totic value. However, side- wall normal stress shows an 
opposite trend. As can be seen from fig. 13, ct„ increases 
with /ip before it eventually saturates to a final value at 
larger y^. Considering the fact that the plic angle of re- 
pose B grows with both [i^ and [iw, one can deduce that 
both wall-grain and grain-grain frictional forces increase 
the stability of the sandpile but in two different ways. 
Increasing [iw while [ip is kept constant, enhances total 
wall frictional forces as the maximum wall friction HwFn 
{Fjy can be each of the wall normal forces Ft, Fg, or Fy) 
might increase even if a^ decreases. On the other hand, 
by increasing grain-grain friction coefficient, the horizon- 
tal normal stresses increase, as suggested by fig. 13. These 
also leads to the strengthening of the maximum frictional 
force experienced at each side walls. 

Of particular interest is the way friction coefficients 
affect the behavior of the Janssen coefficient. We have 
chosen four different values for /iu, : 0.05, 0.3, 0.5 and 1.0, 
and calculated the variation of k with /Zp for each case. 
The results are summarized in a unique diagram in fig. 14. 
From this figure we find that putting aside the case fj.p = 
0.05, the value of k starts to grow rather rapidly at small 
values of ^p, to a final maximal value which is itself an 
increasing function of y^i in all curves. Our simulations 
also show that for the case /Xp = 0.05 the value of k is 
almost insensitive to the variation of wall-particle friction 
coefficient. This might be assigned to the isostaticity of 
the granular packing when ^p is too small, similar to the 
phenomenon observed in Ref.7. Besides, the saturation 
of K with fip does not show a simple trend: while the flat 
part of the first three curves happens to start at /.ip > 
fiw , it begins at a value of /Zp much smaller than fi^j for 
the fourth curve (/Zp = 1.). This observation resembles 
to the results of Ref.7 too, which demonstrate the force 
ambiguity does not follow a monotonic variation with ^p. 

It should be noted that for a more comprehensive study 
of the effect of friction coefficients on k we need to extend 
the range of variation of ^p and fi^ on the logarithmic 
scale. But the problem is that for small values of k a heap 
hardly formed. In fact, a more suitable system for such 
experiments is a granular pack formed in a silo instead of 
a granular heap. Because the formation of a well defined 
pack is always guaranteed even if the frictional forces are 
not significant. 
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FIG. 12: Variation of bottom pressure ab witli grain number 
A'^ in a cell of thickness w = 6 for different values of ^p-. 
fip = 0.1 (□); fip = 0.3 (O); ^p = 0.5 (A); fip = 0.7 (V); 
/ij, = 0.9 (+). 
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FIG. 13: Variation of horizontal stress a^ with grain number 
A'' in a cell of thickness w = 6 for different values of ^p-. 
lip = 0.1 (□); fip = 0.3 (O); ^p = 0.5 (A); /ip = 0.7 (V); 
Mp = 0.9(+). 



V. APPROACHING TO 3-D SANDPILES 

The effect of the cell thickness w on the stability of 
a confined sandpile has been the subject of both exper- 
imental and numerical studies^ril. Careful experiments 
and measurements reveal that the presence of front and 
rear walls increases greatly the angle of repose 6 as well as 
the maximal angle of movement 9m (the maximal angle 
that a granular medium can reach when carefully tilted). 
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FIG. 14: The dependence of Janssen coefficient k on grain- 
grain friction coefficient Hp for different values of /in,: /in, = 
0.05 (D); fi^ = 0.3 (O); fi^ = 0.5 (A);^^., = 1.0 (+). 



These observations show that by increasing the cell thick- 
ness, the angle of repose, 0{w), decreases and it eventu- 
ally settles down to an asymptotic value 0oo when the 
cell thickness becomes very large compared to the size of 
grains. The observed phenomenon is universal and does 
not depend on the formation history of the pile. 

To investigate the effect of w on the wall normal 
stresses, wc conducted a series of simulations on cells 
with different thicknesses, where all the other system pa- 
rameters ad boundary conditions are as listed in Tablel. 
The results for at and cr„ are presented in figl5 and fig. 16. 
These curves indicate again that normal stresses are go- 
ing to saturate with mass for all the cell thickness w but 
with different rates. An interesting feature is that al- 
though (Th grows with cell thickness at smaller values of 
w, but the grows stops when w becomes large compared 
to a single grain size. Meanwhile, the behavior of normal 
stresses at lateral walls is different. In fact the values 
of ay and as (not shown) decrease with w at smaller cell 
thickness but again it seems both of them become almost 
If independent when w is large enough. These behaviors 
arc very similar to the variation of 9 with w^dlr2A. In 
both cases there seems to be a transition from quasi two 
dimensional systems when the cell is just few grain diam- 
eters across, to three dimensional case where the thick- 
ness of the container is much greater than a single grain 
size. 

Let us consider now the behavior of k, the the frac- 
tion of bottom stress a^ which is redirected to each of 
the side walls, with the cell thickness grows to higher 
values. In fig. 17 we have presented our calculation of 
K and 9 for cells with different w, all of them contain 
40000 grains. The other parameters are the basic values 
listed in Tabid. As the figure shows, for small values of 
If the major fraction of at is converted to side-wall nor- 
mal stresses making the contribution of frictional forces 
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FIG. 15: The variation of bottom pressure ab with grain 
number A'^ for various cell thickness w: w — 4(0); w = 6(0); 
w = 8(A); w = 10(V); w = 12(*); w = 14(+); w = 18(*) 



FIG. 16: The variation of front-rear stress, av with grain 
number TV for various cell thickness w: w = 4(n); w = 6(0); 
w = 8(A); w = 10(V); w = 12(*); w = 14(+); w = 18(*) 



more pronounced at these small scales. However when 
the spacing between front and rear walls become large, 
compared with the size of constituent grains, k and con- 
sequently the normal stress exerted by front or rear wall 
approaches constant values. In fact, the behavior of n is 
very similar to the way 6 varies with w. This is a new 
concept and should be considered in physical modelling 
of wall effects on the angle of repose of confined heaps ( 
as well as angle of movement), especially when the spac- 
ing between front and rear walls is of the order of few 
grain diameters. 



VI. CONCLUSIONS 

We have used 3-D MD simulations to investigate the 
effects of a container's walls on the stability of sandpiles 
formed by pouring dried spherical grains into a container 
with small separation between front and rear walls. We 
observed the Janssen's effect, namely the saturation of 
global pressure at the bottom of the container, although 
the geometry is different from the traditional silos where 
the standard Janssen analysis is applied. The same be- 
havior has been observed for front or rear wall's normal 
stresses and the vertical slab supporting the pile. Increas- 
ing the size of the container we witnessed the saturation 
effects again but at larger number of grains. We saw 
that for larger cell thickness, the saturation of normal 
stresses happened at larger values of mass, with stress at 
the bottom plate of the container increases with w and 
the opposite trend for ay. 

We have also studied the dependence of the wall nor- 
mal stresses with wall friction coefficient /i^, and with 
grain friction coefficient fip and found that increasing 
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FIG. 17: Dependence of Janssen coefficient k (squares) and 
angle of repose (filled squares)on the cell thickness w with 
fj.p = fiw = 0.5. 



both of these parameters causes a reduction in normal 
stress felt at the cell's bottom. Meanwhile, we saw that 
the mechanisms of the observed reduction in apparent 
weight while a friction coefficient is increasing, depends 
on which of ^p and fiw is subject to change. While 
the results of simulations show that by increasing wall- 
grain friction coefficient ^y,, the horizontal stress a^ ap- 
proaches smaller values, the changes arc such that the 
overall value of ^wFn which is the maximum of wall 
friction becomes larger and as such a stable heap with 
greater 9 forms. On the other hand, when /x^ is fixed 
and /ip grows, we saw an increase in the calculated val- 
ues of a„, which means the wall friction limit ii^Fisi, is 
again growing with jip. 



We have estimated the Janssen's coefficient k and 
found that not only the grain properties, but also the 
cell geometry and wall- grain interaction affect this quan- 
tity. Our observations suggest that kappa is affected by 
cell thickness w when w is less than about 8d and there- 
fore, the conventional assumption of k being only depen- 
dent to the granular material microscopic characteristics 
is under question at small w. We believe that this is 
a new concept and should be considered in theoretical 
treatment of granular material at small scales. The vari- 
ation of K, with jjLp when other experimental parameters 



are kept constant is also interesting. Our observations 
indicate that k grows almost linearly with ^p for small 
values but reaches to an asymptotic value at larger values 
of /Xp. Meanwhile, when ^p is not too small, increasing 
/i-ii, enhances the pressure transmission from bottom to 
the vertical walls, making k, an increasing function of /iu,. 
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